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Abstract. We discuss the Casimir effect between two identical, parallel slabs, emphasizing 
the role of dissipation and temperature. Starting from quite general assumptions, we analyze 
the behavior of the Casimir entropy in the limit T — > and link it to the behavior of the 
slab's reflection coefficients at low frequencies. We also derive a formula in terms of a sum 
over modes, valid for dissipative slabs that can be interpreted in terms of a damped quantum 
oscillator. 



1. Introduction 

In recent years, the Casimir force |[T] has gained the status of a mesoscopic force, attracting 
considerable attention from physicists and engineers who work with nano- or micro- 
electromechanical devices [21 |3] ID. In particular in view of these applications, many 
theoretical efforts are addressing more realistic configurations. But also for fundamental 
applications, a precise knowledge of experimentally relevant imperfections is valuable to 
tighten upper limits ||5]|6|. The Casimir effect originally describes the force between two 
parallel, perfectly reflecting slabs placed in vacuum at zero temperature, and it is quite clear 
that the generalization to a realistic experimental setup roughly articulates in three steps: 

(i) realistic materials, (ii) non-zero temperature, and (iii) other geometries. Step (i) was 
essentially achieved by Lifshitz in 1955 |7|, and within the same theory, addressing step 

(ii) was apparently trivial. Surprisingly, up to now, the role of dissipation in the temperature 
correction to the force between two dissipative metallic slabs is unclear and at the center of 
hottest controversy around the Casimir effect (ll |9] [TOl [III El El El El- Lifshitz adopts 
a "matter" point of view to calculate the Casimir effect building on a stochastic (Langevin) 
version of the Maxwell equations where quantum fluctuating cuiTents are responsible for the 
fields. Although the inclusion of dissipation is quite natural in this theory, it becomes rather 
complicated when we adopt a "field" point of view. The problem is that the eigenfrequencies 
of the field, uj„i, become complex with absorbing mirrors, so that the expression used by 
Casimir 



is meaningless. We use the symbol [• • -lo^ to denote the renormalization of the sum over 
mode branches to, subtracting the limiting value of the square bracket when the miiTor 
distance L ^ oo. Casimir's expression ([T]i is directly generalizable to partially reflecting 
(dispersive) mirrors Iil6i , but it reveals to be rather tricky in the dissipative case. Among 
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others, Langbein IfTTl [T8l has suggested that one should take the real part of the complex 
energy, but his heuristic arguments remain unsatisfactory. 

In this paper we discuss two different aspects of including dissipation into the Casimir 
effect. We start with a rather general viewpoint on the Casimir entropy at finite temperature, 
recover previous results and give a criterion to classify them. We also show how Eq.([T]i must 
be modified when complex frequencies enter the theory, in order to take consistently into 
account the dissipation phenomena. 



2. The Casimir entropy 

The current controversy arises when calculating the correction to the Casimir force at finite 
(nonzero) temperature. Although the calculation was implemented in Lifshitz's pioneering 
paper and the problem was already recognized by Schwinger et al. [19] , it has been carefully 
examined only recently. An orthodox application of Lifshitz's formula for metal has lead to 
some doubts concerning the way one describes their optical properties. Indeed, using the 
Drude model to describe the optical response of a metallic mirror leads to a completely 
different behavior of the TE and TM polorizations, in stark contrast to perfect reflectors 
(Dirichlet boundary conditions) that, by definition, perfectly reflect both polarizations. This 
phenomenon is fairly striking at high temperatures, but it leads to significant thermal 
corrections also, in particular for the Casimir entropy, in the low temperature limit. In fact, 
one important thermodynamical principle (the Nernst heat theorem) states that when we deal 
with an equilibrium situation the system entropy must go to zero in the limit T ^ 0. 

The final result of Lifshitz's paper is an expression for flie Casimir free energy as a sum 
over Matsubara frequencies: 

j_ oo 

^=^E E'-.9n--'-) (2) 

p,k n— 

where p labels the polarization (TM, TE) of an electromagnetic wave with wavevector k 
parallel to the slab surface. The primed sum has a prefactor 1/2 in the first term, as usual. The 
function g^{uj, r) = In [l — r)^e^^''^] and we use [t is the first Matsubara frequency] 

= fc2 - t^2/c^ T = 2nkBT/h (3) 

Note that we have allowed for a dependence on temperature of the reflection coefficients 
{oj, t). The Casimir entropy can be derived from the previous expression by differentiation. 
Because the series converges uniformly in the parameter t we can interchange the sum and 
the derivative. In the following, we drop the sum over the quantum numbers p and k and set 

kn = 1: 



2TTdT 



5' = -^^ = - ^'[g{inT,T)+mTgu,{mT,T) +Tgr{inT,T)] (4) 

n=0 

where we have defined 

guj{uJ,T) = dujg{uj,T) and ^^(c^, r) ee r) (5) 

Of course the first two terms in (|4|l are more significant than the third, which is identically zero 
when the reflection coefficient does not depend on temperature. We ignore the third term for 
the moment and suppress the second argument for simplicity. The sum over n is re-arranged 
as 



^inr 

n=0 



giinr) - g(i{n ~ 1)t) 
g„(mr) : 



2-3(0) (6) 



Casimir energy and entropy between dissipative mirrors 



3 



Inside the brackets, we use the Taylor expansion of g to get 

g(\{n - l)r) - giinr) + iTg^(inT) = ^ g^k {itit) 

Therefore the entropy per mode can be rewritten as 

5* = - 



fe>2 



fc! 



E, , (-it)'' 1 , , 
ng^u{inT)—^^ :^5(0) 



, n>0, k>2 



(7) 



(8) 



Using the analytical properties of the function g, we can express the fc'th derivative in w by a 
contour integral along a path C in the upper half complex plane (the path must circle around 
each z = inr) and get 




g{z)n{-iTf dz 1 
c — (z-inTf+^ 27ri 2^^ ' 

^ n>0,fc>2 ' 



E 



(9) 



Exploiting the Euler-MacLaurin formula fSol we have {Bm are the Bernoulli numbers] 



fe-2 



inr) 



fc+i 



.,k-l 



1 



fc(fc - 1) 



B2m — 

m>l 



2m- 1)! /T 



(2m-l)!fc! \z 
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E- 

n>0 

When this is substituted into Eq.(|9|l, we get an exact series expansion of the Casimir entropy 
in terms of the temperature and of the derivatives of g{z) si z = Q. 

Let us consider now the limit r — > 0, still ignoring the explicit temperature dependence 
of the reflection coefficients. The sum over fc in Eq.® is approximated by its first term. We 
then find that the entropy goes to zero: 

'1 / 9{z) dz 1 
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This result is the same as the one obtained by Bostrom and Sernelius [S] [21] when 
the dissipation rate in the Drude model does not depend on temperature. The present 
demonstration has the advantage to apply to a broad class of reflection coefficients that do 
not depend on temperature, of course. The first terms of the expansion in t are given by 



S 
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where we express the derivatives in terms of the imaginary frequency ^ = —ioj. In the case 
of the Drude model, fig[T] shows that for the TM polarization, g{i£,) monotonously increases 
with ^. Hence, g^{0) is close to zero or positive. The TE polarization behaves differently: it 
starts from a value g(0) close to zero and shows first a negative slope before increasing and 
going to zero as ^ — > 00; this implies a negative contribution to the entropy. This behavior is 
due to the vanishing of the reflection coefficient r™(w 0) for any finite k. In the general 
case, the sign of the low-temperature entropy depends of course on the balance between the 
derivatives of g for both polarizations; it is therefore model-dependent. 

Nevertheless, we can make some general considerations, which are valid for a large class 
of metals at least as long as the local description for the reflection coefficients is applicable. 
At high frequencies, as long as they are below the plasma frequency, both TM- and TE- 
reflection coefficients are close to unity because the dielectric function is large, giving a 
good approximation to the perfect conductor limit e = 00. This behaviour continues down 
to zero frequencies in TM-polarization because here, electric fields are dominant and are 




Figure 1. Behavior of g{i^) as function of imaginary frequency lo = for the polarizations 
TE and TM. The Dmde dielectric function is used. See the text for more detail. 



efficiently screened by the (Dmde) metal. Only for frequencies much larger than plasma 
frequency ujp, the reflection coefficient goes quickly to zero. Hence, we expect a small and 
positive entropy contribution from the TM mode at low temperatures. In TE-polarization, the 
ultraviolet behavior is the same, but at low frequencies the field is prevalently magnetic and 
penetrates into the material causing that the reflection coefficient rapidly decreases to zero 
II22II . This strongly suggests that the TE contribution to the Casimir entropy is negative and 
larger in magnitude than its TM counterpart. Our analysis also explains why we have always 
a positive entropy for the plasma (non dissipative) model: the UV limit is unchanged, the 
behavior of both reflection coefficients is qualitatively the same, and both do not vanish at 
a; — > 0. A metal described through the plasma model thus behaves qualitatively similar to a 
superconductor at sufficiently small frequencies. We can conclude that a negative value of the 
Casimir entropy is a general feature of the Lifshitz formula when the TE reflection coefficient 
goes quickly and monotonously to zero in the limit of zero frequency. 

Let us consider now a situation where the reflection coefficients r{uj, r) depend on 
temperature. In this case it could happen that the limits r ^ and ^ do not commute. 
More specifically, 

lim r(ar, r) limr(0,T) lim g(aT, t) ^ lim g(0,r) (13) 

where a is some non-zero complex constant. This means that g{oj, r) is no longer an analytic 
function in the whole upper half complex plane when w and r are not independent and we 
approach zero along the path uj = ar. To illustrate this with more detail, let us consider the 
Drude model and the Fresnel form of the reflection coefficients: 

e{uj,T) = l- (14) 

where j{t) is the (temperature-dependent) dissipation rate. The quantity that plays a crucial 
role is uj'^e{uj, r) for which it is easy to see that 

\im{aT)^e{aT,T) ^ lim lim ti;^e(ci;, r) (15) 

if 

lim7(T)/T = (16) 

This happens for the example considered in Ref.|23| where at low temperatures, j{t) oc t^. 
In this particular case, limT-^o r{aT, r) — rpiasma(O). As an immediate consequence, we can. 
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in the general case, no longer perform the previous manipulations for all n. In fact, we have 
to treat n = separately. 

The problem can be handled defining a new function g which is identical to g for all 
nonzero Matsubara frequencies, but does not feature the discontinuity discussed above. In 
the case of the Drude model this can be easily done with the substitution 5(0, t) t-» g{Q) = 
limT^o g{aT, t). The function g so defined has the same analytical properties as a function 
g with no temperature dependence, and we can therefore repeat the previous analysis. The 
substitution has to be made without quantitatively changing our expressions. This means that 
we have to rewrite the expression for the entropy given in Eq.© as follows 



5* = - 



^'[giiTiT, t) +inTg^{mT,T) +Tgr{mT,T)] 



n=0 



(17) 



For the first two terms in the sum we can repeat the previous analysis and show that they 
vanish in the limit t ^ 0. The last term in the sum can be approximated by 



, dg{inT,T) 



n=0 



n=0 



97 



(18) 



which is equal to zero for an at least quadratic temperature dependence of 7(t). In this case 
the low-temperature limit of the Casimir entropy is 

m-m ^ 

2 



< 



(19) 



where we have restored the sum over the wavevectors. Inspection shows that this is nothing 
but the result of Bezerra et al.li23J. 



3. A sum over modes in the dissipative case 

One of the difficulties in generahzing Eq.([T]i to a dissipative system is the concept of a "mode" 
because this represents a state of constant energy which can be well defined only for a closed 
system where dissipation does not occur The naively constructed complex mode frequencies 
actually do not represent the real modes of our system. In fact, in quantum mechanics a 
dissipative system is by definition an open one f24l l25ll26ll27ll28ll29l . i.e. a system coupled a 
macroscopic system, generally a thermal bath. The bath contains a continuum of degrees 
of freedom into which the system energy can be irreversibly dissipated. A formula like 
Eq.([TJ is applicable only if we consider the modes of the total system — which is, of course, 
cumbersome or plainly impossible in practice 130]. It is natural to ask, however, whether 
it may be useful to apply the famihar picture of a damped oscillator in terms of a complex 
frequency, which is just the pole of the Green function for this system. Here we are going 
to show that this compact description can be applied to the Casimir energy as well, with the 
complex poles going over into conventional modes in the non-dissipative case. 



3.1. Sum over poles 

Our starting point will be the Lifshitz formula for the Casimir force between two slabs with 
finite thickness |l7l|3ll[32 33| 

lm[lnDliu)]l (20) 



F — — 9l / did coth 
4-Tr ^ — ' 
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where the frequency integral is just above the real axis (Im w > 0) ensuring retarded response 
functions. The environment is at the scaled temperature r (see Eq.([3]l). We recall the 
renormalization, [• ■ J^, with respect to a pair of infinitely distant mirrors. In the following, 
we drop again the sum over p and k. Finally, the relevant dispersion function is 



Ol In = K 



1 



~2kL 



1 - r{uof 



-2kL ' 



(21) 



with r{uj) being the reflection coefficient of the mirrors (supposed identical for simplicity). 
We used the reality of the optical response, r{uj)* = r(— w*), to extend the integral ( l20b to 
the negative frequency axis. After a few manipulations, the Casimir energy can be written 



E 



2 ^ 



duj 
2^ 



uj coth 



Im[9^1ni?(a>)]^ (22) 



The integrand of this equation goes to zero sufficiently fast for \u)\ — > oo so that we can close 
the integral in the complex plane along a path which encloses the lower half-plane in the 
clockwise sense.. 

The function D is constructed in such a way that its zeros are the complex "mode" 
frequencies we are looking for. But, D also shows branch cuts along the real axis that 
prevent a straightforward deformation of the integration contour These branch cuts come 
in two types: the first corresponds to the mode continuum incident from outside the cavity; 
the second is connected with the dissipative phenomenon and arises from the continuum of 
modes describing the bath. Following Ref. 1 16 1, we eliminate the branch cut of the first type by 
placing our system into a perfectly reflection box with a length a > L. The second type of cut 
is avoided by analytically continuing D from the upper half to the lower half of the complex 
frequency plane. It is exactly because of this step that the zeros of D, otherwise located on a 
unphysical sheet, become apparent |30|. With these changes, D turns into a function Da that 
shows, instead of a branch cut, a set of discrete zeros that adds to the ones of D. One can 
proceed in the usual way with Da and take the limit a ^ oo at the end of the calculation. 

All functions involved in the integral (|22] | are now meromorphic functions |[34l and can 
be expanded over their singularities as follows 



coth 



t/tt 



n— — o< 



d^lllDaiiu)=J2 



1 



(23) 



The zeros of Da must satisfy Imuim < because of causality and passivity. Moreover, 
the dispersion function inherits the reality condition, Da{uj)* ~ Da{~uj*), of the reflection 
coefficients. Therefore, the zeros occur in pairs, ujm and — cl'*„, on both sides of the imaginary 
axis. 

With a little algebra and exploiting the residues theorem one can show 



ujfji coth 



Re 



J L- 



EE- 

_ m n—0 



2nL0^ 



(nr)2 + 



where we have already taken the limit a 

L 

Im " 



^ oo. Note also the following identity 

duj uj Im [di^ In D]^ = 



(24) 



(25) 



than can be easily shown because of parity. This sum rule ensures that in the sum over n in 
Eq.([24]i. the summands go like l/n^ for n ^ oo. 
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As a check let us consider the non-dissipative case: the frequencies are then real, the 
second term in Eq.(l24li vanishes, and the Casimir energy can be written as 



E = 



2 ^ 



LUm {2nT{uJ.m) + 1) 



(26) 



where ^^(wm) is just the mean number of thermal photons in the frequency mode Wm- In the 
Umit of zero temperature this number goes to zero and we get the expressions of Ref. ll35ll36l . 

In the dissipative case, we focus on the zero temperature limit (r 0): the hyperbolic 
cotangent in Eq.(l24]i goes to unity, while the second term becomes an integral so that we have 







n^Re 



d£_ Rc 



2^ 



The integral converges at the upper limit thanks to the sum rule dZST l and we can write 

L 



E= - 



2 ^ 



Re 



2ia;,i 



In- 



A 



(27) 



(28) 



where the meaning of the constant frequency A will become clear soon. This is the main 
result of this section: we have found the Casimir energy between dissipative mirrors, at zero 
temperature, as a sum over cavity mode frequencies. The expression is valid for a generic 
(causal) dielectric function and for mirrors with arbitrary thickness. One sees immediately 
that the non-dissipative case reduces to the standard sum over real mode frequencies. 

To interpret the two terms in Eq.(l28Tl. observe that for each collection of quantum 
numbers p, k, m, we deal with the average energy of a quantum damped harmonic oscillator, 
as discussed in detail by Nagaev and Biittiker (yT\. The first term of the energy is the familiar 
zero point energy, renormalized by the coupling to the environment (see also Ref.|29|). The 
second part of the energy is the most interesting one. It arises because the ground state of the 
oscillator is no longer an eigenstate of the whole oscillator+bath system. (This can be traced 
back to antiresonant terms in the coupling Hamiltonian.) As a consequence, the oscillator 
energy is subject to quantum fluctuations, whose average leads to an extra energy with a 
logarithmic behavior. Within this analogy, the constant A takes the role of a cutoff frequency 
for the coupling to the bath. In fact, for a single oscillator, the cutoff is needed to renormalize 
some quantities like the momentum fluctuations ll37l . It is interesting to note that for the 
Casimir energy, the sum rule dZSl l removes the cutoff dependence in the final result. 



4. Conclusions 



We have explored the impact of dissipation and of temperature on the Casimir effect. On 
fairly general grounds, we have demonstrated that the Casimir entropy vanishes at T ^ 
(the Nernst theorem is satisfied) when the material properties (reflection coefficients) behave 
continuously as both frequency and temperature are decreased to zero. A discontinuity in the 
plane of these variables leads to a non-vanishing, typically negative entropy. We have argued 
that this can be traced to the behavior of the TE reflection coefficient at low frequencies. 
Both results are in agreement with previous work addressing specific models. The complex 
mode frequencies that occur between dissipative slabs are finally a useful tool to compute 
the Casimir energy as a generalized sum over "modes", although the mode concept has to 
be treated with care in this intrinsically open system. Our expression Eq.(l28]l shows that the 
Lifshitz formula (that we started with) does not reduce to the "recipe" of taking the real part 
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of a complex sum over modes, as stated in Refs. IfTTl [TSl . The latter sum deeply changes 
its structure if dissipation is allowed for, while the Lifshitz expression remains formally the 
same. We have argued that the logarithmic correction that we find has a clear meaning from a 
physical point of view, by drawing the analogy to a damped quantum oscillator As it happens, 
the result is independent of the cutoff frequency for the coupling to the environment. 
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